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We present a general relativistic framework for studying gravitational effects in quantum 
mechanical phenomena. We concentrate our attention on the case of ultra-relativistic, spin- 



?— j ' 1/2 particles propagating in Kerr spacetime. The two-component Weyl equation with general 

^2 ! relativistic corrections is obtained in the case of a slowly rotating, weak gravitational field. 

, Our approach is also applied to neutrino oscillations in the presence of a gravitational field. 

The relative phase of two difi'erent mass eigenstates is calculated in radial propagation, and 
1/^ ' the result is compared with those of the previous works. 

I \ 

CN ' SI. Introduction 

> ■ 

I It had been thought that physical phenomena in which gravitational effects and 

O ■ quantum effects appear simultaneously are far beyond our reach, before Colella, 

, Overhauser and WerneJD conducted an elegant experiment using a neutron interfer- 

\^ I ometer. (This kind of experiment is called a COW experiment.) The COW experi- 

■ ment was the first experiment that measures the Newtonian gravitational effect on 

O , a wave function. This effect and its detectability were first suggested by Overhauser 

and Colella,!^ and the next year the effect was verified by Colella et al.^ Although 
' their analysis, which was based on inserting the Newtonian gravitational potential 

. . \ into the Schrodinger equation, was very simple, this experiment is conceptually very 

. ^ ' important in the history of quantum theory. 

^ ■ Recently, gravitational effects on another physical phenomenon, neutrino os- 

^ \ dilations, have been much discussed. '^''^''^''^'''^'•Si The COW experiment and this 

■ ■ ' phenomenon have the common aspect that gravitational effects appear in the quan- 

tum interference. However, there are some differences between the two. In the former 
case, the spatial spread of the wave function plays a significant role, whereas in the 
latter case, the existence of different mass eigenstates and linear superposition are 
important. Another important difference is that the particle is non-relativistic in 
the former case, whereas it is ultra-relativistic in the latter case. 

It seems that a controversy concerning gravitationally induced neutrino oscilla- 
tion phases has arisen. Ahluwalia and Burgard^' state that the phases amount to 
roughly 20% of the kinematic counterparts in the vicinity of a neutron star. Nev- 
ertheless, the definition of neutrino energy and the derivation of the phases are 
not clear in their original paper.^ On the other hand, other groupP^''^'^^''^ have 
obtained similar results for radially propagating neutrinos (the results seem to be 
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different from tliat in Ref. i^). However, the authors of Ref. 4J assume that differ- 
ent mass eigenstates are produced at different times. This assumption seems to be 
of questionable vahdity because the relative phase between the two different mass 
eigenstates initially becomes arbitrary. All of these papers, except Ref. (71), are based 
on a previous work,'^ in which the classical action is taken as a quantum phase. 
Therefore, effects arising from the spin of the particle are not considered in these 
papers. On the other hand, the authors of Ref. use the covariant Dirac equation, 
but they also calculate the classical action along the particle trajectory in the end. 

We provide another framework, different from those of the previous works, for 
studying general relativistic gravitational effects on spin- 1/2 particles with non- 
vanishing mass, such as massive neutrinos. (The experimental confirmation demon- 
strating that neutrinos have nonzero mass has not yet been obtained. However, 
recent experimental reportfl^ seem to suggest neutrinos to be massive.) We do not 
merely calculate the classical action along the particle trajectory, but start from 
the covariant Dirac equation. Our approach allows us to discuss the effects of the 
coupling between the spin and the gravitational field. In particular, we consider the 
propagation of the particle in the Kerr geometry, by which the external field of a 
rotating star can be described. We perform our calculations in a slowly rotating, 
weak gravitational field, and derive the two-component Weyl equations with correc- 
tions arising from the non-vanishing mass and the gravitational field. Furthermore, 
we discuss neutrino oscillations in the presence of the gravitational field. 

The organization of this paper is as follows. In ^ we assume that the external 
field of a rotating object is described by the Kerr metric and discuss the covariant 
Dirac equation in this field. In ^ we derive the Weyl equations with general rel- 
ativistic corrections for an ultra-relativistic particle. The application to neutrino 
oscillations in the presence of the gravitational field is discussed in ^ Finally, we 
give a summary and conclusion in ^ 

§2. Covariant Dirac equation in Kerr geometry 

In this section, we consider the covariant Dirac equation in the presence of a 
gravitational field arising from a rotating object. We derive an equation for the time 
evolution of spinors, which describe particles with spin-1/2, in the last part of this 
section. 

2.1. Cavariant Dirac equation 

To begin, we briefly review the covariant Dirac equation .^"^"^ 

The natural generalization of the Dirac equation into curved space-time gives 

lp- = 0, (2-1) 

where are the covariant Dirac matrices connected with space-time through the 
relations 



d 



r. 



mc 



7V + 7V = 25^^ 



(2-2) 
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and is the spin connection. The spin connection is determined by the condition 

^ - r^^ix - r,iu + i.r, = o. (2-3) 

We now introduce the constant Dirac matrices 7*^"^ defined by 

^(«) = e(«),7^ (2-4) 
where e*^"^ is the orthogonal tetrad satisfying the relation 

g^.u = riabe^^le^'l (2-5) 

{r]ab = diag (c^, —1, —1, — l))- Using these constant Dirac matrices, the spin connec- 
tion is expressed as 



1 



9ux e JV^e.,^ (2-6) 



where the square brackets denote the usual commutator.*) 
2.2. Space-time 

Here, we discuss the gravitational field arising from a rotating object. We now 
assume that the external field of the rotating object is described by the Kerr metric. 
If we restrict ourselves to a slowly rotating, weak gravitational field up to the first 

order in the angular velocity, which is related to the Kerr parameter a, and the 
Newtonian gravitational potential (p = —GM/r, respectively, the line element is 
given by 



ds^ ~ { 1 + 2 ) c'^dt'^ + (xdy - ydx) dt 



- (1 - 2 ^) {dx^ + dy^ + dz^) , (2-7) 

where a is expressed in terms of the mass M and the angular momentum J of the 
gravitational source: 

Assuming that the rotating object is a sphere of radius R with uniform density, we 
have 

where uj denotes the angular velocity of this object. (If the rotating object deviates 
from a sphere, or has an inhomogeneous density distribution, then the numerical 
factor 2/5 might be changed by a factor of order unity.) 



*^ We have ignored a term proportional to the unit matrix. 
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2.3. Equation for time evolution of spinors 

Next, we turn our attention to time evolution of spinors. The covariant Dirac 
equation (|2-H) has beautiful space-time symmetry. However, in order to investigate 
the time evolution of spinors, we must break the symmetrical form of this equation. 
For this purpose, we now use the (3+1) formalism. In the (3+1) formalism, the 
metric ga/3 is split as 

goo = N^ -7ijN'N\ (2.10a) 
goi = -li,N^ = -Ni, (2-lOb) 

gij = -lij, (2-iOc) 

where is the lapse function, A^* the shift vector, and ^ij the spatial metric on the 
3D hypersurface. Furthermore, we define 7*-' as the inverse matrix of jij. Using the 
metric ()2-7|) derived in the last section, we can write the lapse function, the shift 
vector and the spatial metric in the following way: 



N = c{l + ^], (2-11) 



^ = ^^^y^ (2-12a) 
AGMR^ 

= (2.12b) 

A^^ = 0, (2-12c) 



7i, = (l-2^)(5i,-. (2-13) 



Furthermore, we choose the tetrad as 



1 /V* 



= (0'V))' (2.15) 
where the spatial triad e^^'^ is defined as 

7ii e(fe) e^i^ = 6ki. (2.16) 
From Eqs. (|2.11j) - (|2.13|) . we derive 

e(o? = l-|, (2-17a) 

^(0) = — ^^"^y^ (2- 17b) 

2 AGMR^ 

Ho) = 5e2^3 (2.17c) 

6(0) =0, (2.17d) 
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^m= 1 + ^ V- (2-18) 



Using our choice of the tetrad, the covariant Dirac matrices 7" are written as 

/ = 7("^e(,? = 7(°)|, (2-19) 

y = 7(")e(4 = -7(o)^Ar^+70-)e(^.5. (2-20) 

Hence the covariant Dirac equation (|2-H) can be written as 
d 

ih—^ = 
dt 

(^j(o)j(^)cNe^ji - N'^ {pi + ihEi) + ihro + -f^^^mc^N ^, (2-21) 

where p^ is the momentum operator in flat space-time. This equation describes the 
time evolution of spinors. If we adopt the Weyl representation as the constant Dirac 
matrices in this equation, then for massless particles in flat space-time we derive the 
well-known Weyl equations 

d 

ih—ip = ±ca ■ pip, (2-22) 
where ip denotes two-component spinors. 

§3. Ultra-relativistic limit 

We now restrict our attention to the ultra-relativistic limit, which means that the 
rest energy of the particle is much smaller than the kinetic energy in the observer's 
frame. In particular, we expand the energy of the particle itself up to O [m?c'^/pc) . 

Here, we obtain the ultra-relativistic Hamiltonian up to the order of interest by 
performing a unitary transformation similar to the Foldy-Wouthuysen-Tani (FWT) 
transformation . ' > 1 

First, following the discussion of Ref. I14|) . we redefine the spinor and the Hamil- 
tonian according to 

^' = H' = j'/^H^-'/\ (3-1) 

where 7 is the determinant of the spatial metric: 

7 = det(7ij-)- (3-2) 
Since the invariant scalar product is 

(1^,^) = ji^^^d^x, (3-3) 

under this redefinition the scalar product comes to assume the same form as in flat 
space-time: 

{'ilj',(p') = [W^'d^x. (3-4) 
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It is sometimes convenient to adopt this definition of tlie scalar product. 

Next, we perform a unitary transformation to derive the ultra-relativistic Hamil- 
tonian which is the "even" operator up to the order of our interest. From this, we 
have 

H' = UH'U^ 

Hr 
Hl 



+ 



O 



or 0\^,u' 



(3-5) 



The Dirac spinor is also divided into each of two-component spinors as 



(3-6) 



where the subscript R and L denote the right-handed and the left-handed compo- 
nents, respectively. 

We consider the left-handed component. We find that the equation for this 
component is given by 



1 + 



, 1 GM ^ ^ 1 
+ p- (pp^ + 2^^L ■ S 



^^'2 3 p 2 



_ cr ■ p 

cp — 
p 



1 fAGMR^ 



V" 



1 + 



5^3 

1 

4? 



• (L + 5) + 



QGMR^ 
5r^ 



5 • [r X (r X a;) 



=2# + ^P- (PP-P- (PP^ 



+ 2 



p 

1 GM 



p- 

L S 



p' 

rpO p 



2 S 

m c cr ■ p 

2p p 



^ 1 2 2 1 
+ -m C -rr 



A 



crp /rp 1 

P^ P^ P'^ 



where 



A 



AGMR^ 



u-{L + S) + 



QGMR^ 



5 • [r X (r X u>)] 



(3-7) 



(3-8) 



5^3 5^5 

The details of the calculations are given in AppendixEI From this, we find how the 
spin-orbit coupling, the coupling between the spin and the rotation of the gravita- 
tional source, or the coupling between the total angular momentum and the rotation 
is coupled to the non- vanishing mass. 

In radial propagation, the orbital angular momentum vanishes. Therefore, in 
this case, only spin effects coupled to the rotation appear. If we set cj = 0, then 
there is no spin effect in radial propagation. This consequence is consistent with the 
results of previous work.^ 



§4. An application 



In this section, we consider an application of the two-component equation derived 
in the last section to neutrino oscillations in the presence of a gravitational field. In 
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neutrino oscillations (see, e.g., R.efs. ITTj) .ITS Il .119 ]) for analysis in flat space-time), the 
most important point is the phase difference of the two different mass eigenstates. 
Hence we now concentrate on the phase shift of the particle. 

4.1. Neutrino oscillation in Kerr space-time 

We derive the phase shift directly from the two-component equation derived 
in the last section. Furthermore, for simplicity, we consider the radial propagation 
(r-direction), in which the spin-orbit coupling vanishes. 

We now regard terms arising from the non- vanishing mass and the gravitational 
field as perturbations. Then Eq. 1)3 -71) for the left-handed component is considered 
to be 

ih^^i^L = {H^L + ^HL)^L. (4-1) 

where Hql denotes the unperturbed Hamiltonian Hql = —ccr ■ p, and AHl the 
corrections arising from the non-vanishing mass and the gravitational field. 
Here we assume that the spinor V'L is given by 

i;L{x,t) = e"^^'^^OL{x,t), (4-2) 

where -i/'OL satisfies the equation 

d 

^h—ljjQL {x, t) = HoL'fpoL {x, t) (4-3) 

ot 

Substituting Eq. into Eq. (|^ and using Eq. we obtain 



1 



^ = --J AHidt. (4-4) 

In order to derive the gravitationally induced phases practically, we assume that 
corresponding to the left-handed component, ^oL satisfies the relation 

^-^i)QL{x,t) = -iljQL{x,t). (4-5) 

p 

Furthermore, we here replace the g-numbers in AHi with c-numbers. This is a kind 
of semi-classical approximation. From this, excluding spin effects, we derive the 
phase 

1 rri^c^^ 



where we have considered the case that the neutrino is produced at a space-time 
point A{tA,rA) and detected at a space-time point B{tB,rB)- We now concentrate 
on the term related to m^, because neutrino oscillations take place as a result of 
the mass square difference. Let the two different mass eigenstates have common 
momentum p and propagate along the same path. Then the relative phase A^Pij of 
the two different mass eigenstates, and {i^j), is given by 
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Next, we discuss the spin-rotation coupling. In a similar way, we replace the 
g-numbers in terms concerning the spin of the particle in Eq. (|3-7() with c-numbers 
again. From this kind of semi-classical approximation, we find that the spin-rotation 
term coupled to m? vanishes. However, the spin-rotation effects in higher order 
terms may survive. This fact implies that there is no influence of the spin-rotation 
coupling on neutrino oscillations, at least up to the order related to m^. 

Consequently, in a radially propagating case, we obtain the phase difference (|4-7j) 
between the two different mass eigenstates as a final result. 

4.2. Comparison with previous works 

Finally, we compare the result obtained above with those of previous works. For 
this purpose, we consider the Schwarzschild limit (i.e., a? — > 0), which leads to the 
result ()4-7() again. 

First, we assume the "background" neutrino trajectory as that of radial null 
geodesies: 

= ds'=(^l + 2^) c'dt^ - (^1 - 2^) dr\ (4-8) 

Then we obtain 

Hence, if we transform the integral ()4-7() with respect to t to that with respect to r, 
the relative phase ^^ij is given by 



The second term in the round brackets corresponds to the gravitational correction as 
indicated by Ahluwalia and Burgard.l^ Indeed, under the assumption that the tetrad 
component of the radial momentum p(^) = &(^^Pr is constant along the trajectory, we 

can obtain the same expression as in Ref . . (Note that Ahluwalia and BurgarcP^ 
assume P(r)C as the energy of the neutrino.) 

Next, let us see whether our result (|4-7() reproduces the other form of the re- 
sults.'^'t^''^ From the mass shell condition g^^p^Pu = rn^c?, p is related with the 
energy E{= ptc) in the following way: 

pc={l- 2^ E+[0 (m^) or O (0^)] . (4-11) 

Under the assumption that E is constant along the trajectory, we finally obtain 



1-3 



221 

2h .1,^ E 



Amf-c^ 

which is the same as the result obtained in previous works.'^'^''^ In this sense, 
our result Eq. (|4-7() contains both of the previous expressions. Our analysis here 
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clearly shows that the controversy arising due to the discrepancy between the re- 
sults of Ahluwalia and Burgard^' and other author^ '^^''^ is simply due to different 
assumptions concerning constancy along the neutrino trajectory. 

§5. Summary and conclusion 

We have studied the general relativistic effects of gravity on spin- 1/2 particles 
with non-vanishing mass. In particular, we have considered particles propagating in 
the Kerr geometry in the slowly rotating, weak field approximation. By performing 
a unitary transformation similar to the FWT transformation, we have obtained the 
two-component Weyl equations with the corrections arising from the non-vanishing 
mass and the gravitational field from the covariant Dirac equation. The Hamilto- 
nian clearly shows how the spin-orbit coupling, the spin-rotation coupling or the 
coupling between the total angular momentum and the rotation is coupled to the 
non- vanishing mass. 

Furthermore, we have discussed an application of the two-component equations 
to neutrino oscillations in the presence of a gravitational field, and we have derived 
the phase difference of the two different mass eigenstates in radial propagation. It 
is worth mentioning that our result contains both of the previous expressions. We 
have shown that the controversy arising from the discrepancy between the results of 
Ahluwalia and Burgard^' and other authorP^'^''^ is simply due to different assump- 
tions regarding constancy along the neutrino trajectory. Moreover, as seen in the 
transformation of the integral variable, the gravitationally induced neutrino oscilla- 
tion phases arise from the modification of the propagating distance. Indeed, we found 
that the gravitational correction term comes out in this variable transformation. 

We have not applied our approach to a non-radially propagating case in detail 
in this paper. However, it is of interest whether the spin-orbit coupling affects the 
neutrino oscillations. This will be the subject of further investigation. 

Although it seems difficult to provide verification of these effects with current 
experimental detectability, we believe that investigation of systems in which both 
quantum effects and gravitational effects come into play is important. Progress in 
technology may make the verification of such effects possible. 
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Appendix A 

Derivation of Ultra- Relativistic Hamiltonian 



A.l. Components of spin connection 

The spin connection is given by Eq. H2-6|) : 

1 r 



gxa e(„f V^e 



(A-1) 



It is convenient to introduce the following 4x4 matrices similar to the Pauli spin 
matrices: 



Pi 



/ 
/ 



P2 



-il 
il 



P3 



/ 
-/ 



(A-2) 
(A-3) 



where I is the 2x2 unit matrix. These matrices satisfy the relations 

PiPj = ^ij + i^ijkPk, 

where Sijk is the Levi-Civita antisymmetric tensor (ei23 = +1)- If we adopt the Weyl 
representation as the constant Dirac matrices, then we have 



(A-4) 



where a are the Pauli spin matrices. 

Using the above quantities, the components of spin connection, up to the order 
of interest, are given by 



i^^o = P'i ■ {P<4 
2c 



+ 



1 



u ■ S -\ — — S • [r X (r X u)] 



ihFi 



ihFo 



ihr?, 



c 

h 

2? 

ih 3GMR^ 



+ 

h 

2? 



P3 



5r^ 5r^ 

[-2xy Gi + (x^ - y^) - yz (T3 



5^6 



-LO 



',3 0"! - 0-3) 



ih SGMR"^ 



+ ^P3- 



5r^ 



-co 



[(x^ - y^) cri + 2xy 0-2 + zx 0-3] 



2c2 



,1 0"2 



\2 CTl) 



ih SGMR'^ , 

+ — P3 — =-5 — [-yz CTl + zx a2\ 
c or 

A. 2. Unitary transformation 

The Hamiltonian defined in Eq. (|2-2ip is given by 
1 



(A-5) 



(A-6) 



(A-7) 



(A- 



H = psca ■p + P3 



-ccr ■ p 



+ 2^ C<T • p 
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5r^ 



u-{L + S) + 



5r^ 



S ■ [r X {r X cj)] 



2 2 V 

+ pimc + pimc -5-. 



(A-9) 



Moreover, the Hamiltonian redefined by Eq. (|,'-{-H) is tlien 



H = psca ■ p + { CO- ■ p ^ + CO- ■ p 



+ 



1 



AGMBf' 
5^3 



u-{L + S) + 



QGMR^ 
5^5 



S -[r X {r X u))] 



2 2^ 
+ pimc + pimc 



(A-IO) 



where terms proportional to pi and p2 are "odd", and those proportional to ps are 
"even" . 

Next, by performing a unitary transformation similar to the FWT transforma- 
tion, let us derive the ultra-relativistic Hamiltonian for the left-handed component. 
Here we divide the unitary transformation into several steps. First, we use the 
unitary operator 



Ui = exp ( ip2-mc 



' 2 c^p 

which is introduced to eliminate the odd term p\m(? . Using the useful formula 



(A-ll) 



e^^He-'^ 



H + i[S,H] + t, [S, [S, H\] + ^ [S, [5, [S, H]]] + 



and the relation (|A-3p . we obtain the transformed Hamiltonian 
UiH'ul 

= p^ca • p + p3 ccr • p — + — c<T • p + -^A + pimc 



(A-12) 



Pi-mc 
1 



C<T • p 



C(T ■ p 



C<T • p 



g2p2 g2 g2 c^p^ 



+ 



CCT • p 



C<T • p 



g2p2 g2 g2 (.2^2 



ccr • p 



2 f C<T ■ p A A CIT ■ p\ 1 24Cf-p 

^ T'/^ C I I j yf-^^ ^t^y~l/ 



1 2 4 

P3-m c 



1 



^2p2 g2 



ca ■ p + ccr ■ p 



1 



g2 g2p2 



ccr • p 



+ 



ccr • p 



g2p2 g2 g2 g2p2 



1 2 4^^ 1 

-m c 



g2 g2p2 



ca ■ p A ca ■ p 1 A 

2 „ n — „ n + 



g2p2 g2 g2p2 



g2p2 g2 



where A is given by Eq. H3-8() : 



QGMR^ r~i I / \^ 
u-{L + S) + — • [r X (r X u)] . 



(A-13) 



(A.14) 



5^3 5^5 

Second, in order to eliminate the second line in Eq. ()A-13|) . we use the unitary 
operator 



U2 = exp 



1 2 ( CCT-p 



-ip2-mc 



+ 



ca • p 



g2p2 ^2 g2 ^2p2 



(A-15) 
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Using this unitary operator, we obtain 



U2UiH'ulul 



psca ■ p + p-i [ CO- ■ p + CO- ■ p j + + pimc 



1 A ccr ■ p 



ip2-mc 



'- 2 4 

P3-m c 



^2p2 g2 g2 g2p2 



ca ■ p + ccr ■ p 



1 2 A^cr p 
1 f ca -p q 



^2p2 g2 



g2 g2p2 



g2p2 g2 



+ 



(p ca ■ p 



1 24/^ 1 
-m c 



g2 g2p2 



ccr • p j4 c<T • p 1 A 

2 „ n — ^ n + 



g2p2 g2 g2p2 



g2p2 g2 



(A-16) 



Finally, we use the two unitary operators f/3 = e*"^^ and f74 = e*"^"*, where ^3 and ^4 
satisfy the relations 



i [S^^p-ica ■p] = - pimc^^, 
i [S'4, p-ica ■p] = - ip2^rnc^ 



ca ■ p A A ca ■ p 

g2|j2 ^2 g2 g2p2 



(A-17) 
(A-IS) 



We here assume the existence of these unitary operators, which make the remain- 
ing odd terms vanish. (We need not find the concrete forms of these unitary op- 
erators, because extra terms arising from these unitary transformations are higher 
order terms.) Using these unitary operators, we obtain the transformed Hamiltonian 
UH'W which is even up to the order of our interest. 



/ + _ / _ </> (A -A 1 

UHU^ = p'ica - p + pslccr-p-^+^ccr-pl + -^A 



1 2 4ccr -P 

2 c^p 
1 

■I 

A 1 



1 2 4 

— P3-m c 



2 4 

-m c 



0=9 -^ccr ■ p + ca ■ p-^ -^r^ 

c'^p c^p 



ca ■ p (p (p ca ■ p 

^2p2 g2 g2 ^2p2 



ca ■ p A ca ■ p 1 A 

2 „ o — „ o + 



^2p2 g2 ^2p2 



ip,p2 g2 



(A.19) 



where U is given by C/ = U4U3U2U1. From this, we can derive Eq. (|3-7() by a simple 
calculation. 
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